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ABSTRACT 

We derive the basic physical parameters of the field double-mode RR Lyrae star 
BS Com from its observed periods and the requirement of consistency between the 
pulsational and evolutionary c onstraints. By using the current solar-scaled horizontal 
branch evolutionary models of iPietrinferni et alj (|2003 ) and our linear non-adiabatic 
purely radiative pulsational models, we get M/M Q = 0.698 ± 0.004, log(L/L ) = 
1.712 ±0.005, T cS = 6840 ± 14K, [Fe/H]=-1.67±0.01, where the errors are standard 
deviations assuming uniform age distribution along the full range of uncertainty in 
age. The last two parameters are in a good agreement with the ones derived from 
the observed BVIq colours and the updated ATLAS9 stellar atmosphere models. We 
get T off = 6842 ± 10K, [Fe/H]=-1.58 ± 0.11, where the errors are purely statistical 
ones. It is remarkable that the derived parameters are nearly independent of stellar 
age at early evolutionary stages. Later stages, corresponding to the evolution toward 
the asymptotic giant branch are most probably excluded because the required high 
temperatures are less likely to satisfy the constraints posed by the colours. We also 
show that our conclusions are only weakly sensitive to nonlinear period shifts predicted 
by current hydrodynamical models. 

Key words: stars: individual: BS Com - stars: variables: RR Lyrae - stars: abun- 
dances - stars: evolution - stars: fundamental parameters 



1 INTRODUCTION 

Double-mode radial pulsators play a very significant role 
in the study of pulsating stars. This is because we can 
determine all fundamental stellar parameters just by us- 
ing a few observed parameters and relying on the rather 
solid theory of linear radial pulsation (assuming of course, 
that the mode identification is correct - a question of- 
ten raised for typically low-amplitude pulsators whereas 
rarely bothered with in the c ase of high /moderate-amplitude 
ones) . In a series of papers l|Kovacs fc W alker 1999; Kovacsl 
2000a. tj) we investigated the consistency of the method 
by using the observed colours and assumed cluster metal- 
licities in deriving distances for globular clusters and for 
the Magellanic Clouds. These distances have been con- 
firmed by independent studies, for example by direct 
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HST parallax meas u remen ts of Ivan Leeuwen et alj (|2007f ) 
and iBenedict et al.1 (120071). by Hipparcos-based 5 Scut i 
PL-relation of McNamara. Clementini fc Marconi (|2007f ) . 
by mixtures o f par allax and interferometric methods of 
iFouque et al.| l|2007l) and by B aade-Wesselink analyses 
I: IKqvIcs1|2003|) . 



(Storm ct al. 2004; Kovacs 2003). In yet another applica- 
tion. iBuchler fc Szabdl 1 2007 ) put limits on the metallicities 
of Magellanic Cloud beat Cepheids through stellar evolution 
constraints. 

Double-mode radial pulsators are also viable for non- 
linear hydrodynamical modelling. Although early attempts 
have already shown that sustained double-mode pulsation 
is possible with purely radiat ive hydrodynamical code s (at 
least for RR Lyrae stars, see lKovacs fc Buchlerl 0*9931 ). the 
modelling was put on a better sounding physical ground only 
after time-depen dent convection was more carefully included 
(for Cepheids by Kollath et al.lll998l. for RR Lyrae sta rs by 



lFeuchtingerlll99S and ISzabo. Kollath fc Buchlerlliooi ) Al- 
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though these models are capable of producing stable double- 
mode pulsation in the appropriate neighbourhood of the 
physical parameter space and the model light curves are in 
reasonable agreement with the observed ones, star-by-star 
modelling and application of more stringent observational 
constraints are still lacking. 

Large-scale surveys aimed at microlensing and vari- 
ability studies have already yielded a considerable increase 
in the number of known double-mode variables. In par- 
ticular, before 2000, only four double-mode RR Lyrae 
(RRd) stars were identified i n t he Milky Wa y (see 
iGarcia-Melendo fc Clemendll997l and lKovacsH2001bl for ad- 
ditional statistics). At the moment of this writing we have 
26 RRd stars known in the Galactic field according to the 
literature, mostly due to the NSVS and ASAS full-sky vari- 
ability surveys, and many more in t he Galactic Bulge, due to 
the OGLE survey (lMizerskil l2003l; iMoskalik fc Porettill2003l 
and IPieulski. Kolaczkowski fc KopackH l2003t ) . 

The field RR Lyrae star BS Com (a = 13 h 34 m 39 s , 8 = 
+24° 16' 38") was originally put on the observing list of the 
60 cm autom ated telescope of the Konkoly Observatory (see 
Sodor 2007 for details on the ongoing photometric survey), 
du e to the suspected Blazhk o behaviour of this star reported 
bv lClementini et alj dl995bl) . While gathering data, the star 
was analysed bv lWilsl (|2006l ). using the NSVS database. He 
found that BS Com was an RRd staifl 

This paper describes the results of the frequency anal- 
ysis of our data and the pulsational/evolutionary analysis. 
The result of this latter, purely theoretical study will be 
confronted with additional constraints posed by the effective 
temperature and metallicity obtained from our multi-colour 
data. We show that the purely theoretical approach yields 
stellar parameters in agreement with the colour data and 
can be used for very accurate parameter estimation once we 
are able to pose an age limit for the object. 



2 PHOTOMETRY 

2.1 Observations and data reduction 

We obtained BV(RI)c photometric observations of BS Com 
with the 60 cm automated telescope of the Konkoly Obser- 
vatory, Budapest. The telescope is equipped with a Wright 
750 x 1100 CCD provid ing a 17" x 2 4' field of view (detailed 
parameters are given bv lBakos|[l999h . The observations were 
carried out between 2005 March 1 and 2006 April 17 in 

3 seasons, on a total number of 35 nights. Observational 
statistics are shown in Table [T] BS Com magnitudes were 
measured relative to BD +24 2598 by aperture photome- 
trjB GSC 00454-00454 and USNO-B1 1143-0206728 served 
as check stars. We used standard packages for data 
reduction procedures. Coefficients for linear transformation 
of the instrumental data into standard magnitudes were de- 
rived each year. 

We note however, that selected the ~ 1 yr alias of 

the fundamental period. 

2 Photometric data are available electronically on-line at CDS. 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 1. BS Com observational statistics 



date 


nights 




data points 




(JD - 2450000) 




B 


V 


Rc 


Ic 


3431 - 3467 


19 


568 


566 


561 


543 


3744 - 3762 


7 





148 


149 


149 


3824 - 3843 


9 





256 


214 


243 


total: 


35 


568 


970 


924 


935 



Table 2. Standard BVIq magnitudes 



Star 


V 


B — V 


V-I C 


BD+24 2598 


10.608 


0.603 


0.662 


BS Com 


12.755 


0.286 


0.433 


errors: 


±0.0015 


±0.0016 


±0.0032 



Notes: Magnitude averages are given. The errors are standard 
deviations of the mean differences between BD +24 2598 and 
HD 117876. 

BS Com was tied to the standard photometric sys- 
tem in B, V and Ic colours through the comparison star 
BD+24 2598. The magnitudes of BD +24 2598 were mea- 
sured relative to HD 117876 (V = 6 m 082, B — V = 0?968 
V - R = ffi°522, V - I = r?017, see iHoltzman fc Nations] 
1984), the closest standard star to BS Com. The observa- 
tions were carried out near two subsequent culminations, 
under very good weather conditions. Because of the bright- 
ness of HD 117876, a large fraction of the incoming light had 
to be blocked in order to avoid saturation. We reduced the 
effective aperture of the telescope to the ~ 7% of its origi- 
nal size by placing a circular diaphragm at the front end of 
the telescope tube. Because of the ~ 24' angular distance 
between HD 117876 and BD+24 2598, it was not possible to 
observe the stars in a common field of view. In order to keep 
our photometric precision, the two distinct fields were ob- 
served repeatedly. We obtained 5, 10, and 9 frames in B, V, 
and Ic, respectively, for each of the two stars. Data were cor- 
rected for atmospheric extinction. The standard magnitudes 
and their formal statistical errors are shown in Table [2] For 
BS Com, magnitude averages are shown, calculated over the 
complete light curve solution given in Sect. 12.21 

2.2 Fourier analysis: the double-mode solution 

We performed a standard frequency analysis on the V 
data by discrete Fourier-transformation. We derived Po = 
0.48791 d and Pi = 0.36307 d for the periods of the funda- 
mental and the first overtone modes, respectively. Besides 
the frequencies of the two pulsation modes, 5 harmonics and 
8 linear combinations of the frequencies corresponding to the 
above periods were identified in the amplitude spectra. In or- 
der to increase the detection probability of potentially exist- 
ing additional componen ts, we applied th e spect rum averag- 
ing method described bv lNagy fc Kovacsl (120061 ). Due to the 
averaging of the frequency spectra in the V(RI)c colours, 
the noise of the residual spectra decreased by a factor of 
~ v3, down to the ~ 2mmag level. Except for the above 
two components, their harmonics, and linear combinations, 
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no additional components wer e found. More details on the 
frequency analysis are given bv lDekanvl(|2007f) . Our 15-term 
Fourier-sum fit gives a full description of the light curves 
within the observational accuracy. The Fourier-parameters 
and their errors are given in Table [ATI 



3 PHOTOMETRIC ABUNDANCE AND 

TEMPERATURE FROM B,V,I C COLOURS 

Knowledge of the metallicity and temperature is very im- 
portant in the determ ination of the physical param eters 
of double-mode stars (fCol Il99ll ; iKovacs et all Il992l and 
lKovacsll2000bl ). There are two former estimates of the metal- 
licity of BS Com. ISmithl (|l990l ) obtained two measurements 
of AS = 6.0 and 8.2. The second value was taken at 
lower brightness than the first one, so this, being closer 
to minimum li ght yields [Fe/H] = —1 .53 on the scale of 
IJurcsikl <| 19951 ). Iciementini et al.l (|l995bl ) acquired one spec- 
trum of BS Com with a resolution of 1.4 A and a signal-to- 
noise ratio of ~ 20. They derived a very low abundance of 
[Fe/H] ~ —2.0 from t he equivalent widt h of the Calcium K- 
line. As explained by Iciementini et al.l (|l995bl) . they took 
the spectrum 'at minimum light', according to their photo- 
metric V time series consisting of 97 points. However, taking 
into consideration the double-mode behaviour of BS Com 
yet unknown at that time, the real pulsational phase of the 
aforementioned observations is ambiguous. (Unfortunately, 
the pulsational parameters of BS Com derived in this pa- 
per do not enable us to predict the light curve accurately 
enough at the time of these observations.) Therefore, we de- 
cided to make an independent estimation of the metallicity 
and also the effective temperature of BS Com based on our 
photometric data. 

The idea of deriving heavy element abundances from 
magnitudes measured in proper photo metric bands is not 
new. In his early study, ISturchl (|l966l ) used UBV magni- 
tudes around minimum light to estimate t he line bla nket- 
ing of RRab stars. In a similar application, iLubl ll 19791 ) em- 
ployed the reddening-free Walraven [B — L] colour index for 
metallicit y determination. Also, we refer to a more recent 
paper by iTwarog. Anthonv-Twarog fc De Lee] (|2003h . who 
obtained more precise abundances for open clusters using 
intermediate-band photometry. 

The method represents the most economical approach 
in terms of getting a first order estimate on [M/H]. In a com- 
pletely formal approach, the method of photometric abun- 
dance determination comes from the following idea. In the 
static stellar atmosphere computations the models are de- 
termined by fixing basic input parameters, such as chemi- 
cal composition, convective parameters (mixing length over 
atmospheric scale height, turbulent velocity, overshooting, 
etc.), gravitational acceleration g and effective temperature 
T e ff. If we fix the admittedly poorly known convective pa- 
rameters and assume that solar-type abundance distribution 
is, in general, a good approximation of the true distribution, 
we are left with three parameters only: overall metal abun- 
dance [M/H], log g, and T e g. Then, it is obvious that once we 
have an independent estimate on log g, then we can compute 
both T e ff and [M/H] from two colour indices, assuming that 
they depend differently on the above two parameters. Be- 
cause of the availability of the average B,V,Ic magnitudes 



of BS Com and because the B — V colour index relativel y 
strongly depends on [M/H] (see e.g. lKovacs fc Walkerll999T l. 
whereas V — Ic is mostly sensitive to T c g , it is useful to ex- 
amine the possibility of determining [M/H] and T c g through 
the B — V and V — Ic colour indices. 

Because of the poor representation of the B — V, 
[M/H], logg, and T c s dependence by a single linear or low- 
order polynomial formula, we decided to use directly the 
grids available for stellar atmosphere models. We employ 
quadratic interpolation in order to get accurate model colour 
indices for each set of physical parameters. More precisely, 
having the value of log g fixed, we scan the parameter range 
in B— V , V — Ic and find the best [M/H], logT c ft parameters 
that minimise the following function: 



V = [logT cff (B - V) - logT off (V - Ic)] 2 
+a B . v [(B-V) ohB -(B-V)] 2 
+a 2 V -i c [(V - / c )ob s - (V - Ic)] 2 , 



(1) 



where T c g(B— V) and T c g(V — Ic) are the T ff values at the 
respective colour indices, the subscript 'obs' means the ob- 
served values, whereas the scanned values are denoted with- 
out subscripts. The weights ob-v and otv-i c are set equal 
to 0.33 and 0.25, respectively, in order to account for the pro- 
portionality of log T aH to 0.33(5- V) and 0.25(V-I C ) in for- 
mula e obtained by linear regressions fsee lKovacs fc Walker! 
1999). The introduction of the last two terms in T> is nec- 
essary, because of the high sensitivity of the method for ob- 
servational noise. 

W e emplo y the stellar atmosphere models of ICastellil 
(1999, see also ICastelli. Gratton fc KuruczHl997T ) computed 
with the opaci ty distribution func t ions b ased on the solar 
abundances of iGrevesse fc Sauval (| 19981 ) (grids quoted as 
'ODFNEW at Kurucz's web siteQ). Th ese new grids dif- 
fer fro m the previous ones as described in IPietrinferni et al.l 
(2004). All models have a microturbulent velocity of 
2kms~ , a convective mixing length parameter of 1.25, 
scaled solar heavy element distribution and no convective 
overshooting. For the computatio n of log g we used a similar 
formula to that in lKovacsl (|2000al ). adapted to the appropri- 
ate parameter regime of RR Lyrae stars: 



log fl = 2.47 - 1.23 log P 



(2) 



We note two possible sources of ambiguity when using 
Eq. (flj) for the determination of T e g and [M/H]. First, the 
zero points of the various colour — T B g calibrations (e.g. in- 
frared flux method, IRFM versus calibrations based on Vega 
or Sun-like stars) may be slightly different. In earlier applica- 
tions w e used zero points defined by IRFM (see e.g. iKovacsl 
2000a). For pulsating stars, there is also a problem of the 
differen ce between the static and pulsation -phase-averaged 
colours (|Bono. Caputo fc Stellingwerfll 19951 ). Since the sizes 
of these effects are still debatable, we opted to use the stel- 
lar atmosphere data bv lCastellil (|l999l ) without applying any 
corrections to the magnitude-averaged colours. As we shall 
see in the subsequent test, the derived metallicities are in a 
good agreement with the spectroscopic ones, indicating that 
the theoretical colours are fairly accurate. 

4 http: / /kurucz. harvard.edu 
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Table 3. Reddening dependence of the photometric [Fe/H] and 
T c ff values of BS Com. 
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Figure 1. [Fe/H] p h values computed from the BVIq colours 
plotted against the spectroscopic [Fe/H] sp values taken from the 
literature for 24 RRab stars. Error bars denote standard devi- 
ations for each star, obtained by adding uncorrelated Gaussian 
random noise with a = 0^005 to the observed B — V, V — Iq 
colour indices. Open circles show the outlying stars X Ari, SS Leo, 
and VY Ser. 



To test the reliability of the method we collected 24 
fundamental mode RR Lyrae (RRab) stars from the litera- 
ture with known periods, B — V and V — Ic average colour 
indices, and spectroscopic [Fe/H] valued). Details on this 
test sample are given in Appendix [B] In Fig. [T] we show the 
derived [Fe/H] p h values versus the [Fe/H] sp spectroscopic 
metallicities of the test stars. By the exclusion of the three 
stars X Ari, SS Leo, and VY Ser with discrepant photometric 
metallicities (open circles in Fig. 1), the standard deviation 
of the fit decreases from 0.26 dex to 0.15 dex. We consider 
the latter as the 1 a error of the method. Estimates of the er- 
rors of the individual data points were computed by adding 
Gaussian random noise of a — 0^005 to the input colour in- 
dices. The error bars in Fig.[T]indicate that the sensitivity of 
the method to observational errors increases with decreasing 
[Fe/H] (see also Kovacs 2008 for more details). 

We estimated the metallicity and effective tempera- 
ture of BS Com by using logg = 2.85 from Eq. © and 
the average B — V and V — Ic values given in Table [2] 
The amount of reddening at the high Galactic latitude of 
BSCom (/ = 19?00, b = +79?86) is expected to be low . 
Indeed, according to ISchlegel. Finkbeiner fc Davisl l|l998T ). 
E(B — V) = 0.014. This gives a maximum value for the real 
reddening of the star, and therefore poses minimum values 
to its intrinsic colours. The estimated metallicities and effec- 
tive temperatures of BS Com with different possible redden- 
ing values are shown in Table [3] Formal statistical errors are 
10 K and 0.11 dex, computed the same way as for the test 
objects, but using the photometric errors listed in Table [2] 



E(B - V) 


[Fe/H] 


Tea 


0.000 


-1.47 


6774 K 


0.005 


-1.51 


6798 K 


0.010 


-1.54 


6823 K 


0.014* 


-1.58 


6842 K 



* ISchlegel et al] l|l998l l 

4 PULSATION MODELS 

We have generated linear non-adiabatic (LNA), fully ra- 
diative stellar models to derive the fundamental physi- 
cal parameters of BS Com. The pu lsation code is basi- 
cally identical to the one describe d by [ Buchlerl (ll99C ), orig- 
inally develo ped by [Stellingweril (Il975h fsee also Castor! 
Il97lf) . As in iKovacs (2006), we used the Rosseland mean 
opacity values taken f rom l lglesias fc Rogers! (|l996h and 
I Alexander fc Fergusonl (|l994f ) interpolated to the appro- 
priate chemical composition. T he heavy element d i stribu - 
tion is solar-scaled, as given by iGrevesse fc Anders! (|l99lf ). 
The stellar envelopes were densely sampled with 500 mass 
shells down to the inner radius and mass values, where 
Q in = (ftin/iW)(M in /M surf ) = 0.05, with R suvl and M surf 
values at the outermost grid point. This zoning resulted in a 
deep envelope down to (6 — 7) ■ 10 6 K and ensured a sufficient 
(better than 10~ 4 d) stability in the periods (as compared 
to shallower envelopes of Qi n = 0.10). 

We covered the metallicity and temperature ranges well 
beyond the ones derived from the BVIc colours in Sect. [3] 
By scanning these ranges and matching the observed peri- 
ods, we obtained unique solutions for the possible physical 
parameters of the star. These solutions are listed in Table [4] 
To be m ore compatible w it h the evolutionary models de- 
rived by iPietrinferni et al.l |2004), with the change of the 
overall metal content Z, we have also changed the hydro- 
gen abundance X. We note that our best estimated [Fe/H] 
and T c ff values are closest to the item at Z = 0.0006 and 
Tcff = 6850 K. It is noticeable, how sensitive the solution is 
to the metal content. For example, by changing [Fe/H] by 
~ 0.2 dex, the mass changes by 0.1 -0.15 Mq. However, a 
200 K change in T e ff causes a change in the mass less than 
0.06 Mq. We also draw attention to the remarkable stability 
of the gravity and density. While the [Fe/H] change men- 
tioned above causes ~ 15 % change in mass, the same figure 
is less than 0.8% and 0.3% in log g and logp, respectively. 
The effect of temperature change is even smaller. It is also 
important to note that all models listed in the table are lin- 
early excited in their fundamental and first overtone modes. 
Therefore, they are viable for sustained double-mode pulsa- 
tions. More detailed tests by direct nonlinear hydrodynamic 
simulations are beyond the scope of this paper (see, however 
Section [SJ. 



COMBINATION OF THE PULSATIONAL 
AND EVOLUTIONARY MODELS 



5 For solar-scaled heavy element distribution we have [M/H] = 
[Fe/H], 



Stellar evolutionary models play a crucial role in under- 
standing basic stellar physics and developing more coherent 
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Table 4. Pulsation models matching the periods of BS Com 





M 


logL 


R 




log 9 


logp 






X = 


7^.47 




— u.uuuo 




6700 


0.639 


1.647 


4.940 




2.856 


-2.127 


6750 


0.646 


1.664 


4.964 




2.857 


-2.128 


6800 


0.656 


1.682 


4.993 




2.858 


-2.130 


6850 


0.663 


1.699 


5.018 




2.858 


-2.131 


6900 


0.673 


1.717 


5.049 




2.860 


-2.133 






V — 
A — 


7^4R 
U. / 0440 


y 






6700 


0.669 


1.661 


5.021 




2.861 


-2.128 


6750 


0.678 


1.679 


5.050 




2.862 


-2.130 


6800 


0.688 


1.697 


5.080 




2.863 


-2.131 


6850 


0.698 


1.715 


5.111 




2.865 


-2.133 


6900 


0.709 


1.733 


5.143 




2.866 


-2.134 






V — 
. \ — 


O 7^'iA 
U. / Oo4 




— U.UUUD 




6700 


0.732 


1.690 


5.188 




2.872 


-2.132 


6750 


0.743 


1.708 


5.221 




2.873 


-2.133 


6800 


0.754 


1.726 


5.253 




2.875 


-2.135 


6850 


0.768 


1.745 


5.291 




2.876 


-2.137 


6900 


0.780 


1.763 


5.324 




2.877 


-2.138 






X = 


0.7530 


z 


= 0.0010 




6700 


0.872 


1.745 


5.530 




2.893 


-2.139 


6750 


0.885 


1.763 


5.563 




2.894 


-2.140 


6800 


0.901 


1.782 


5.603 




2.896 


-2.142 


6850 


0.914 


1.800 


5.637 




2.897 


-2.143 


6900 


0.929 


1.818 


5.672 




2.898 


-2.145 



Notes: T c ff is given in [K], M, logL and R are measured in solar 
units, g and p are in [CGS]. Observed periods match within 
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models. For example, the 'Cepheid mass discrepancy' prob- 
lem before 1992 appeared as a basic confrontation between 
the 'standard' evolutionary and pulsational theories. A des- 
perate search for the clue of this discrepancy led Norman 
Simon to a plea for th e re- visitation of the heavy element 
opacities llSimon 19821). that has res ulted fin ally in a suc- 
cess m 1992 (|Rogers fc Iglesiaslll992l . see also ISeaton et al.l 
1 19941 ). In another application of the evolutionary results, in 
the current search for extrasolar planets, a generally used 
method for determining the physical parameters of the host 
star is matching stellar isochrones in the p s tar - Tcff plane to 
the observed values (assuming that T e ff and [Fe/H] are re- 
liably determined from e.g. spectroscopy, see lSozzetti et ail 
l2007t ). 

From a completely formal point of view, the method of 
parameter determination by the combination of evolution- 
ary and pulsational models is equivalent to the following 
problem. At a certain pair of periods, the LNA pulsational 
model grid (Table 2} gives stellar mass and luminosity as a 
function of temper ature and metallic i ty. M oreover, HB evo- 
lutionary tracks of IPietrinferni et alj (|2004h also give stellar 
mass and luminosity for arbitrary temperature, metallicity, 
and age values. Therefore, we have the following two vector 
functions sampled in discrete points: 

E: ([Fe/H], Teg , t) — * (M,L) 

P : ([Fe/H], Tcff , Po,Pi) -> (M, L) , 




Figure 2. Stellar mass and luminosity as functions of the metal- 
licity. Age and temperature are fixed at log(t[yr]) = 7.424 and 
T B ff = 6851 K. Continuous and dashed lines show the quadratic 
interpolants in Z for the evolutionary (squares) and pulsational 
(circles) models. For better visibility, the intersections are con- 
nected by a vertical grey line. 



where E stands for the evolutionary models, P denotes the 
pulsational models, and t is the time elapsed from the start 
of the core helium burning at the zero age horizontal branch 
(ZAHB). If we assume consistency between the two kinds of 
models, then the subspace of stellar parameters that fulfil 
the constraints posed by both models is given by the inter- 
section of the two vector functions (E = P). This will be 
(at the most) a ID subspace, i.e. a set of solutions that has 
age t as a free parameter. The [Fe/H] and T c ff ranges corre- 
sponding to our multitude of solutions will then be compared 
with the ones obtained from the photometry in Section [3] 
The method is implemented as follows. 

(a) We construct is ochrones from the s olar-s caled horizon- 
tal branch models of IPietrinferni et al.l (|2004l ). Physical in- 
gredients and details of the model building are described in 
their paper and references therein. Here we note only some 
relevant features. All models are canonical in the sense that 
they have been calculated without atomic diffusion and con- 
vective overshooting. They include up-to-date input physics, 
most importantly, they use recent and accurate results for 
plasma-neutrino processes and electron conduction opacities 
which lead to more accurate He core masses at evolutionary 
stages after the core helium flash. And last, but not least, 
their HB models have a fine mass spacing, particularly suit- 
able for investigating stellar parameters in the instability 
strip. 

(b) From the pulsational solution (PS) we get 'near 
straight lines' for any fixed T e ff in the Z — > (M, log L) 
planes (hereafter ML planes). 

(c) For a given age we can search for the evolutionary so- 
lution (ES) by using the above fixed T e g. These solutions 
yield another set of 'near straight lines' in the ML planes. 

(d) The ES and PS lines may or may not intersect in the 
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Table 5. Physical parameters of BS Com derived in this paper 



[Fe/H] 


Tcff [K] 


M/Mq log(L/L Q ) 


-1.67 ±0.01 


6840 ± 14 


0.698 ±0.004 1.712 ±0.005 


-1.58 ±0.11 


6842 ± 10 





Note: Errors for the ES/PS solution (first line) are the standard 
deviations of the various age-dependent solutions in the full age 
interval of 60 Myr. Items derived from the BV Iq colours (second 
line) are taken from Table [3] errors are formal statistical ones. 



appropriate Z range in the ML planes. If they do, then 
we have two Z values at which intersection occurs. Obvi- 
ously, the consistent ES = PS solution should yield the same 
Z. Since this is generally not the case, the best solution is 
searched for by minimising the difference between the two 
metallicities: 



AZ = \ log Z M -logZ L \ 



(3) 



where Zm an d Zl are the metallicities obtained in the M and 
L planes, respectively. To illustrate the solution obtained at 
a given age and temperature, we show the PS and ES curves 
and their intersections in Fig. 2. It is seen that these curves 
are indeed near linear. Therefore, we can employ quadratic 
interpolation to compute the positions of the intersections. 

(e) By changing the age, with the above method we can 
map the solution as a function of the evolutionary stage. As 
it follows from the above construction, at each stage we have 
different stellar parameters that yield the same periods and 
satisfy evolutionary constraints. The grey-scale plot of AZ 
on the log(t) — T e g plane is shown in Fig. 3. It is remarkable 
how much the solution is limited in a narrow T e s regime. It 
is also observable that we get solutions (i.e. AZ < 10~ 4 ) at 
nearly all ages up to log(t) = 7.77 (« 60 Myr). We could, 
in principle, fix the age by using very accurately (and in- 
dependently) determined T e g. However, the precision of the 
currently available methods for temperature estimation is 
insufficient for posing a strong constraint on the age. 

The ridge line in Fig. 3 selects the best fitting models. 
The age dependencies of the corresponding log L, M, and 
[Fe/H] are shown in Fig. 4. The small range of stellar pa- 
rameters allowed by the solution is remarkable. The stronger 
topological change occurs at logt = 7.6 (~ 40 Myr) just 
after the blueward loop where evolution starts toward the 
asymptotic giant branch. Considerably larger ages seem to 
be less probable because the associated high temperatures 
can be excluded by the ones derived from the colours. From 
the above result we can estimate the parameters of BS Com. 
Assuming equal probability of ages in the range plotted, we 
can compute the distributions of the various physical pa- 
rameters. Although these distributions are not symmetric 
around their most probable values, for the present case this 
effect is of secondary importance and we compute the stan- 
dard deviations around the averages as a first approximation 
of the errors introduced by the effect of age uncertainty. Ta- 
ble [5] displays the so-derived stellar parameters, which are 
compared with the independently estimated T e g and [Fe/H] 
values from Table [3] It is reassuring that there is a comfort- 
able overlap between the two sets of parameters, indicating 




Figure 3. Grey-scale plot of the metallicity difference AZ [see 
Eq. i3t ] on the age — T of j plane. 
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Figure 4. Age dependence of the derived stellar parameters for 
BS Com. 



a good level of consistency among the very different ingre- 
dients of the above parameter determination. 

The employment of the LNA pulsation models in the 
above solution presumes that the theoretical periods are 
applicable to the observed ones that are known to be in- 
herently nonlinear due to the longevity of the observable 
double-mode state. Since our method strongly relies on this 
assumption, it is important to examine if the currently avail- 
able nonlinear models predict some considerable change in 
the asymptotic nonlinear periods with respect to the LNA 
periods. We note that this problem is also present in any as- 
teroseismological investigation if the multimode state is the 
result of some nonlinear phenomenon. However, nonlinear 
modelling is usually not easy to attempt since most of the 
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multimode pulsators are of non-radial, therefore it would re- 
quire full (i.e. 3D) hydrodynamics. We can pose the above 
question only in the case of RRd stars and double-mode 
Cepheids, since these variables are thought to be pure ra- 
dial pulsators and we have successful nonlinear modelling 
available on these variables. 

Data on the nonlinear/linear period difference for 
asymptotic hydrodynamical models are nea rly non-existent. 
In an early assessment l|Kovacsl l2001br i we used the 
somewhat limited and preliminary result obtained by 
iKollath fc Buchlerl l|200ll) . The nonlinear periods were 
longer than the linear ones in the order of ~ 0.001 days. 
The amounts of the period shifts were such that they re- 
sulted in a decrease of less than 0.002 in the period ratio. 
An ongoing, more extended study by Szabo et al. (to be pub- 
lished) basically confirms the above result. As it is detailed 
in Appendix [C] we generated another set of LNA models 
by artificially decreasing the observed periods of BS Com by 
the amount predicted from the models of Szabo et al. After 
repeating the same computation as for the LNA periods, we 
ended up with stellar parameters that were remarkably close 
to the ones obtained with the LNA assumption. The solution 
curves (age vs. stellar parameters) exhibit small overall shifts 
that result in differences between the parameters obtained 
with the LNA and modified LNA solutions. The parameter 
shifts are approximately as follows: +0.01, +0.006, +10 K, 
-0.08, for M/Mq, log(L/L ), T cB and [Fe/H], respectively, 
in the sense of modified LNA minus non-modified LNA. We 
conclude that considering the period shifts predicted by the 
currently available nonlinear simulations does not effect sig- 
nificantly the agreement with the independently estimated 
T cft and [Fe/H]. 



6 CONCLUSIONS 

We performed a pulsational/evolutionary analysis on one 
of the most extensively observed Galactic double-mode 
RR Lyrae (RRd) star BS Com. Our earlier analysis incor- 
porating all of our measurements in the BV(RI)c colours 
has shown that the light variation does not contain any 
other components except for those of the fundamental and 
first o vertone modes together with their linear combina- 
tions (|Dekanvll2007h . Therefore, the presence of possible 
non-radial components is excluded above the ~ mmag level. 
(This may not be the case for AQ Leo, where the analysis 
of the high precision data from the MOST satellite suggests 
the presence of two p ossibly non-radial compon ents at the 
sub- mmag level - see iGruberbauer et all 120071 .) Based on 
the result of the frequency analysis, we decided to derive 
the physical parameters of BS Com using the assumption of 
radial double-mode pulsation. 

Earlier studies on double-mode variables were mostly 
based on the period ratio (Petersen) diagram, which is an ef- 
fective (but not complete) representation of the main (mass 
and heavy element) dependences of the periods (see e.g. 
|Popielski'et~aHl2000h . It is also possible to employ other 
available pieces of i nformation, su ch as colour or overall clus- 
ter metallicity (e.g. lKovacsll2000ah . In this paper we followed 
a purely theoretical approach, where we derived all basic stel- 
lar parameters by matching the evolutionary and pulsational 
properties. In the solutions age remains a free parameter. 



It is important to confront the above theoretical solu- 
tion with other independent pieces of information. For this 
goal, we utilised the observed BVIc colours in deriving 
metal abundance and effective temperature with th e help of 
updat ed solar-scaled stellar atmosphere models by ICastellil 
(1999). We tested the reliability of the method through the 
comparison of the published spectroscopic abundances of 21 
Galactic field fundamental mode RR Lyrae (RRab) stars. 
The standard deviation of the residuals between the spec- 
troscopic and photometric metallicities is 0.15 dex, which 
we consider as a good support for the reliability of the pho- 
tometric method based on the BVIc colours. By apply- 
ing this method on BSCom, we get [Fe/H]= -1.58 ± 0.11, 
Tcft = 6842 ± 10 K, where the errors are purely statistical 
ones. 

Taking into consideration the above independently es- 
timated parameters, we can most probably exclude the late 
evolutionary stages after the core helium exhaustion. This 
is because the high temperature required by the solution is 
incompatible with the temperature computed from the ob- 
served colours. Therefore, in the acceptable set of solutions 
the age spans only an approximately 60 Myr range from the 
zero age horizontal branch till slightly after the start of the 
final redward evolution. Because we do not have independent 
information about the evolutionary status, all of these solu- 
tions are viable. However, the ranges covered by the param- 
eters in the possible time interval are remarkably narrow. 
Furthermore, assuming uniform age distribution, most of the 
values are clumped in an even narrower regime at the ear- 
lier evolutionary stages. By considering this age ambiguity, 
we get M/Mq = 0.698 ± 0.004, log(L/i ) = 1.712 ± 0.005, 
T cff = 6840 ± 14K, [Fe/H]=-1.67 ± 0.01 where the errors 
are 1 a ranges due to the above age interval. 

It is interesting to compare the above values with the 
compilation of the average physical parameters of the RRd 
sta rs found in vario us stellar systems. In Table 2 of the paper 
by iKovacs! l|2001al ) we summarised the result for four sys- 
tems, where the parameter determinations were performed 
without the use of evolutionary models. There are IC 4499 
and LMC that have average metallicities similar to that of 
BSCom. Omitting LMC, due to its very wide population 
range, we see that the average RRd parameters in IC 4499 
(T cff = 6760 K, logL/L = 1.708, M/Mq = 0.755) are re- 
markably similar to those of BS Com. The largest difference 
occurs in the mass, but this is due to the slightly longer 
period of BS Com than the average value of the RRd stars 
in IC4499. Since the other stellar parameters are nearly the 
same, this period difference results in a larger average mass 
for the cluster stars. 

Another indirect check of the consistency of our the- 
oretical parameters of BS Com is the computation of the 
ze ro point of the period-l uminosity-colour (PLC) relation 
of IKovacs fc Walker! (|200 IT ) . Although this relation was de- 
rived from fundamental mode RR Lyrae stars, because of 
its apparent generality, we use it with the period of the 
fundamental mode of BS Com. By using Mbol.o = 4.75, a 
bolometric correction interpola ted to the app ropriate stellar 
parameters of BC = —0.15 (se e Castclli 1999) and a redden- 
ing value given bv lSchlegel et al.l "( 19981 . see Table 3), we get 
— 1.02 for the PLC ze ro point . In c omparing this value with 
the one published by iKovac 1 (|2003l). we see that the above 
value is larger by 0.06, yielding lower distance modulus for 
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the LMC by the same amount (we rec all that t he LM C dis- 
tance modulus with the zero point of lKovacsl (|2003h comes 
out to be 18.55). We think that the above result confirms the 
consistency of the derived physical parameters, especially if 
we add the ambiguity due to the yet to be shown fit of the 
RRd variables to the PLC relation derived from RRab stars. 

We have also tested the effect of a possible period shift 
due to convection and nonlinearity. Based on the extensive 
set of models by Szabo et al. (in preparation) , our conclusion 
is that the effect in the derived parameters is relatively small 
(~ 1.5%). Once this correction is proved to be important 
by future studies, taking it into consideration is fairly easy, 
since the overall period shift seems to be almost independent 
of the particular model parameters. Therefore, there is no 
need to run specific nonlinear models in each case but it is 
quite enough to fit the relatively inexpensive linear models to 
the modified periods based on the above general corrections. 

The remarkably small ambiguity of the derived stellar 
parameters from the purely theoretical method employed 
in this paper brings up several intriguing questions. One of 
these is the compatibility of the theoretical abundance and 
temperature and the ones obtained by spectroscopic and 
various semi-empirical methods (e.g. by the widely used in- 
frared flux method). Accuracies such as 1% (or better) are 
quoted fo r the temperature in the cur rent papers on this 
topic fe.g. iMasana. Jordi fe R ibas 2006). Unfortunately, for 
the chemical composition we still lack this level of accuracy. 
Most of the abundances available for RR Lyrae stars are 
based on some calibration of low-dispersion spectra by a 
low number of high-dispersion ones. Double-mode stars in 
globular clusters are obvious targets for the above-type of 
studies, due to the expected chemical homogeneity and same 
distance. Last, but not least, we note the need of putting 
more stringent constraints on the nonlinear hydrodynami- 
cal models. Our current understanding is that double-mode 
behaviour is caused by the non-resonant interaction of two 
normal modes with a delicate level of physical dissipation 
(understood as the result of convection). Although the global 
observational and theoretical properties are in reasonable 
agreement, the accuracy of the presently available stellar 
parameters is insufficient to lend further support to these 
models. Considering the general astrophysical importance 
of convection and, in particular, the role of that in double- 
mode pulsation, it is clear that pinning down the physical 
parameters is of great significance. 
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APPENDIX A: FOURIER DECOMPOSITION 
OF BS COMAE 

Table IA1I shows the frequencies, amplitudes and phases 
of the 15-term Fourier sum fitted to the light curves of 
BS Com. Phases refer to sine-term decomposition with an 
epoch t = 2453431.0. The Fourier-sums fit the B,V,R c ,Ic 
data with 0.0116, 0.0097, 0.0091, 0.0098 mag rms scatters, 
respectively. These values are in the appropriate proxim- 
ity of the formal errors of the individual photometric data 
points as derived from the IRAF package. To determine the 
errors of the fitted amplitudes and phases, we performed 
Monte Carlo simulations of the light curves. Synthetic data 
were generated by adding independent Gaussian random 
noise to the Fourier solution, with a made equal to the rms 
scatter of the residual data. Errors in Table |A"T1 are the stan- 
dard deviations of the Fourier parameters obtained from 100 
independent realisations of the light curve in each colour. 



APPENDIX B: THE SAMPLE FOR TESTING 
THE BVIc [Fe/H],T eff RELATION 

Here we present the parameters of the 24 RRab used in 
Sect. [3] to test the method of metallicity and temperature 
estimation from the BVIc colours and fundamental peri- 
ods. Table IB1I shows fundamental periods, colour indices, 
interstellar reddenings and spectroscopic [Fe/H] sp values 
taken fro m the literature. W e relie d mainly on the compi- 
lati ons by Kovacs & Jurcsik (1997) for colour indices and 
on lJurcsik fe Kovacsl l|l996h for [Fe/H], except for three 
stars. For the abundances of U U Cet and V440 Sgr we re- 
fer to IClementini et all l|l995al ) and for RVPhe to iJoned 
i| 19731 ). Reddenings were taken from iBlancol (|l992l ). except 
for BB Pup for which we used the E(B — V) value given 



Table CI. Pulsation models matching the modified periods of 
BSCom 





M 


logL 


R 




log9 


logp 






X = 


0.7547 


Z 


= 0.0003 




6700 


0.672 


1.662 


5.026 




2.862 


-2.128 


6750 


0.679 


1.679 


5.050 




2.863 


-2.129 


6800 


0.689 


1.697 


5.080 




2.864 


—2.131 


6850 


0.697 


1.714 


5.105 




2.865 


-2.132 


6900 


0.707 


1.732 


5.137 




2.866 


—2.134 






X = 


0.7546 


Z 


= 0.0004 




6700 


0.705 


1.677 


5.114 




2.868 


-2.130 


6750 


0.715 


1.695 


5.144 




2.869 


-2.131 


6800 


0.724 


1.713 


5.172 




2.870 


-2.133 


6850 


0.733 


1.730 


5.200 




2.871 


-2.134 


6900 


0.745 


1.748 


5.232 




2.872 


-2.135 






X = 


0.7534 


Z 


= 0.0006 




6700 


0.772 


1.706 


5.287 




2.879 


-2.133 


6750 


0.786 


1.725 


5.325 




2.881 


-2.135 


6800 


0.798 


1.743 


5.356 




2.882 


-2.136 


6850 


0.810 


1.761 


5.389 




2.883 


-2.138 


6900 


0.820 


1.778 


5.416 




2.884 


-2.139 






X = 


0.7530 


Z 


= 0.0010 




6700 


0.919 


1.761 


5.633 




2.900 


-2.140 


6750 


0.930 


1.778 


5.660 




2.901 


-2.141 


6800 


0.944 


1.796 


5.694 




2.902 


-2.142 


6850 


0.959 


1.814 


5.728 




2.903 


-2.144 


6900 


0.971 


1.831 


5.757 




2.904 


-2.145 



Notes: T e g is given in [K], M, logL and R are measured in 
solar units, g and p are in [CGS]. Observed periods match 
within 10~ 5 d. The modified periods are P = 0.486448 and 
Pi = 0.362340 days. 

by ISchlegel et all l|l99Sl) . The derived [Fe/H] ph photomet- 
ric abundances for each star are listed in the last column 
of Table |Bl] For XAri, SSLeo, and VYSer, we get dis- 
crepant [Fe/H]ph values. These might be due to system- 
atic errors in their colour indices, which have a particu- 
larly large effect on the estimated [Fe/H] p h at lower metal- 
licities (see Kovacs 2008). We also note that SSLeo and 
VY Ser h ave also been met as outli e rs in various earlier stud- 
ies (e.g. IKovacs fc Kanbuil 1 19981 ; IKovacs fc Walker! 1200 ll ; 
iKovacsl [20031 ) . 



APPENDIX C: EFFECT OF NONLINEARITY 

In Szabo et al. (in preparation) we examine the period shifts 
between the asymptotic nonlinear period of the convective 
models and those of the purely radiative linear non-adiabatic 
models (i.e. the ones we use in this paper). The results are 
based on many different models, including various combi- 
nations of physical and model construction parameters (no- 
tably depth of the inner boundary and number of mass shells 
constituting the stellar model). All models show that both 
modes have longer periods in the asymptotic double-mode 
regime than those in the LNA approximation. The amount 
of period increase is smaller than 0.001 days. The shifts lead 



10 /. Dekdny et al. 



Table Al. Fourier decompositions of the B, V, Rq, Ic light curves of BS Com. 



Identification 


Frequency 


A B 


A v 




A Ic 


fo 


- 


2 04955 


1 64 2fi + 77 


1 97 f)^ -4- 47 


1 ni 1 S + 45 

iwi.io _i_ vj.'-ttj 


7q fiq + fl c;i 


2/0 


— 


4 DQQ1 fl 


26 51 -1-0 70 

iiU.tJJ- _1_ KJ . 1 KJ 


on 87 + 47 


1 7 26 + 46 

J. 1 . £i\J _1_ VJ.'-tU 


1 4 88 + 48 


3/0 




6.14865 


4.43 ±0.78 


1.82 ± 0.47 


1.44 ± 0.45 


2.48 ± 0.47 




fl 


2.75432 


259.43 ±0.76 


201.66 ± 0.46 


161.36 ± 0.42 


123.06 ±0.51 


- 


2/1 


5 50864 


35 55 + 83 


29 08 + 49 


22 76 + 50 

1 'J 1 \ f . ' f\ f 


17 42 + 47 


- 


3/1 


8.26296 


15.57 ±0.68 


12.48 ± 0.46 


9.87 ±0.49 


6.86 ±0.44 


- 


4/1 


11.01728 


3.75 ±0.65 


3.21 ± 0.42 


2.94 ±0.45 


2.48 ± 0.47 


fo 


fl 


4.80387 


62.50 ±0.65 


50.23 ± 0.48 


40.87 ±0.54 


32.82 ±0.54 


fo 


-fl 


0.70477 


36.19 ±0.68 


29.65 ± 0.43 


24.98 ± 0.49 


18.10 ±0.53 


2/0 


fl 


6.85342 


11.59 ±0.78 


9.75 ± 0.51 


8.54 ±0.42 


6.86 ±0.52 


fo 


2/1 


7.55819 


21.67 ±0.74 


20.31 ± 0.47 


15.30 ±0.49 


11.38 ±0.57 


2/0 


2/1 


9.60774 


10.85 ±0.77 


6.08 ± 0.46 


4.99 ±0.41 


3.14 ±0.45 


fo 


3/1 


10.31251 


8.04 ± 0.79 


6.46 ± 0.46 


4.28 ±0.42 


3.61 ±0.50 


fo 


-3/1 


6.21341 


5.09 ±0.83 


2.63 ± 0.50 


3.52 ±0.47 


2.24 ±0.50 


fo 


4/1 


13.06683 


3.35 ±0.56 


2.56 ± 0.45 


1.78 ±0.43 


1.46 ±0.43 
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- 


fl 


2 754320 


q 1 21 + n 004 


q nQfi + n 004 

O.UOU _1_ u.uu^ 


q riQA 4- n n05 


2 ofin + 006 




2/1 




9 occr: _|_ n non 
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3/1 


8.262960 


1.023 ±0.176 


1.562 ± 0.273 


2.202 ± 0.355 


1.433 ±0.260 




4/1 


11.017280 


2.216 ±0.003 


2.200 ± 0.003 


2.172 ± 0.003 


2.104 ±0.004 


fo 




2.049550 


1.560 ±0.020 


1.534 ±0.016 


1.502 ± 0.019 


1.569 ±0.030 


2/0 




4.099100 


0.438 ± 0.046 


0.433 ± 0.035 


0.426 ± 0.048 


0.457 ±0.074 


3/0 




6.148650 


5.553 ±0.180 


5.992 ±0.133 


5.746 ±0.137 


5.727 ±0.169 


fo 


fl 


4.803870 


1.564 ±0.014 


1.576 ±0.008 


1.602 ± 0.013 


1.542 ±0.014 


fo 


-fl 


0.704770 


4.576 ±0.022 


4.616 ±0.015 


4.657 ±0.018 


4.682 ± 0.024 


2/0 


fl 


6.853420 


0.875 ± 0.054 


0.829 ± 0.046 


0.931 ± 0.059 


0.793 ±0.066 


fo 


2/1 


7.558190 


0.234 ±0.035 


0.251 ±0.023 


0.250 ± 0.027 


0.316 ±0.045 


2/0 


2/1 


9.607740 


0.319 ±0.071 


0.159 ±0.077 


0.113 ± 0.071 


0.178 ±0.135 


fo 


3/1 


10.312510 


0.698 ±0.097 


0.359 ± 0.082 


0.102 ± 0.109 


6.266 ±0.114 


fo 


-3/1 


6.213410 


4.988 ±0.162 


5.180 ±0.199 


5.342 ± 0.145 


5.418 ±0.236 


fo 


4/1 


13.066830 


5.980 ±0.164 


5.303 ±0.170 


5.341 ± 0.232 


4.824 ±0.358 



Notes: The light curves are represented in the from of Aq + A\ ■ sin(27r/(t — to) + €>i) + . . . 
with to = 2453431.0. Frequency is given in d , amplitudes and phases arc in mmag and radians, 
respectively. We show Icr formal statistical errors. 



to period ratio values that are lower by up to 0.0015 than pensating effect of changes in T e ff and [Fe/H] (see Sect. 5 
the LNA ones. For the overall relative period shifts we get: for additional details). 



Plna(FU)/P nl (FU) = 0.997 (CI) 

Plna(FO)/P nl (FO) = 0.998 . (C2) 

Assuming that the observed periods of BS Com correspond 

to the nonlinear asymptotic values, the predicted modified 
LNA values are as follows: 

P (LNA, mod.) = 0.486448 (C3) 

Pi (LNA, mod.) = 0.362340 (C4) 

Pi /Po (LNA, mod.) = 0.744869. (C5) 



With these periods we performed the same type of LNA 
survey as described in Sect. 5 for the observed periods. The 
result is shown in Table CI. By comparing the items at 
the same chemical composition and temperature with those 
matching the observed periods, we see that M and log L are 
shifted to larger values by ~ 0.04 and ~ 0.015, respectively. 
However, in the final solution the difference between the pa- 
rameters obtained by the predicted and directly computed 
LNA periods will be considerably smaller, due to the com- 



Pulsational and evolutionary analysis of BS Com 



Table Bl. Properties of the test RRab stars 



Star P FU B-V V - I c E(B — V) [Fe/H] sp [Fe/H] ph 



SWAnd 


0.4422659 





.434 





.544 


0.09 


-0 


.06 


-0.22 


WYAnt 


0.5743365 





.384 


0. 


.549 


0.05 


-1 


.39 


-1.13 


XAri 


0.6511597 





.490 





.724 


0.16 


-2 


.10 


-2.50 


RRCet 


0.5530288 





.371 





.537 


0.05 


-1 


.29 


-1.28 


UUCet 


0.6053409 





.390 





.567 


0.01 


-1 


.38 


-1.33 


WCrt 


0.4120139 





.370 





.466 


0.09 


-0 


.45 


-0.22 


DXDel 


0.4726182 





.454 





.571 


0.09 


-0 


.32 


-0.26 


SU Dra 


0.6604200 





.348 





.519 


0.02 


-1 


.56 


-1.50 


SW Dra 


0.5696710 





.369 





.520 


0.01 


-0 


.95 


-0.93 


RXEri 


0.5872475 





.415 





.585 


0.10 


-1 


.07 


-1.13 


RR Gem 


0.3973082 





.406 





.511 


0.11 


-0 


.11 


-0.23 


RR Leo 


0.4523926 





.334 





.486 


0.03 


-1 


.30 


-1.23 


SS Leo 


0.6263438 





.350 





.488 


0.02 


-1 


.56 


-0.82 


TT Lyn 


0.5974332 





.381 





.550 


0.05 


-1 


.50 


-1.26 


V445 Oph 


0.3970232 





.620 





.772 


0.29 


+0. 


.01 


-0.18 


AVPeg 


0.3903760 





.424 





.528 


0.10 


+0 


.08 


-0.17 


ARPer 


0.4255489 





.676 





.853 


0.35 


-0 


.14 


-0.31 


RVPhe 


0.5964182 





.370 





.535 


0.03 


-1 


.50 


-1.19 


BB Pup 


0.4805468 





.459 





.582 


0.11 


-0 


.35 


-0.47 


V440 Sgr 


0.4774788 





.404 





.568 


0.09 


-1 


.11 


-1.09 


VYSer 


0.7140962 





.377 





.579 


0.02 


-1 


.71 


-2.33 


WTuc 


0.6422370 





.329 





.490 


0.01 


-1 


.37 


-1.39 


TUUMa 


0.5576570 





.364 





.522 


0.05 


-1 


.15 


-1.15 


UUVir 


0.4756062 





.349 





.475 


0.02 


-0 


.60 


-0.66 



Note: Periods are in [days] , colour indices and reddenings are in [mag] . See Appendix [B] for references on 
the sources of data. 



